A PROOF OF THE GLOBAL ATTRACTOR CONJECTURE IN THE 
SINGLE LINKAGE CLASS CASE 
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Abstract. This paper is concerned with the dynamical properties of deterministically modeled 
chemical reaction systems. Specifically, this paper provides a proof of the Global Attractor Conjecture 
in the setting where the underlying reaction diagram consists of a single linkage class, or connected 
component. The conjecture dates back to the early 1970s and is the most well known and important 
open problem in the field of chemical reaction network theory. The resolution of the conjecture has 
important biological and mathematical implications in both the deterministic and stochastic settings. 
One of our main analytical tools, which is introduced here, will be a method for partitioning the 
relevant monomials of the dynamical system along sequences of trajectory points into classes with 
comparable growths. We will use this method to conclude that if a trajectory converges to the 
boundary, then a whole family of Lyapunov functions decrease along the trajectory. This will allow 
us to overcome the fact that the usual Lyapunov functions of chemical reaction network theory are 
bounded on the boundary of the positive orthant, which has been the technical sticking point to a 
proof of the Global Attractor Conjecture in the past. 
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1. Introduction. This paper is concerned with the qualitative behavior of de- 
terministically modeled chemical reaction systems with mass-action kinetics. We will 
provide multiple results pertaining to weakly reversible reaction systems that will al- 
low us to conclude that the Global Attractor Conjecture, the most well known and 
important open problem in the field of chemical reaction network theory, holds in the 
single linkage class case. That is, the conjecture holds when the underlying reaction 
diagram consists of a single connected component. 

1.1. Background and statement of the problem. Natural questions about 
the qualitative behavior of deterministically modeled chemical reaction systems in- 
clude the existence of positive equilibria (fixed points), stability properties of equilib- 
ria, and the non-extinction, or persistence, of species, which are the constituents of 
the system. As the exact values of key system parameters, termed rate constants and 
which we will denote by Kfc, are usually difficult to find experimentally and, hence, 
are oftentimes unknown, it would be best to answer these questions independently of 
the values of these parameters. Building off the work of Fritz Horn, Roy Jackson, and 
Martin Feinberg [HI [H [HI [211 121I23] the mathematical theory termed "Chemical 
Reaction Network Theory" has been developed over the previous thirty-five years to 
answer these types of questions. 

Early work by Feinberg, Horn, and Jackson showed that if a reaction network with 
deterministic mass-action kinetics admits a so called "complex-balanced" equilibrium, 
see Definition 12.81 then there exists a unique complex-balanced equilibrium within 
the interior of each positive compatibility class, or invariant manifold |161 I2H 123) 
(throughout, we will often refer to the invariant manifolds of the systems of interest 
as compatibility classes, or stoichiometric compatibility classes, to stay in line with 
the terminology of chemical reaction network theory). Horn and Jackson also proved 
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the existence of a strict, entropy type, Lyapunov function that gives local asymptotic 
stability of each such equilibrium relative to its compatibility class. Later, Horn, 
Jackson, and Feinberg proved what is best known as the Deficiency Zero Theorem: 
that regardless of the choice of parameters Kk-, a reaction network with deterministic 
mass-action kinetics that is both weakly reversible and has a deficiency of zero must 
admit a complex-balanced equilibrium |13 ( I14 | [T6]. Here, a reaction network is weakly 
reversible if each of the connected components of its reaction diagram is strongly 
connected, see Definition 12.41 and the deficiency of a network is defined in Definition 
12.101 Collecting ideas shows that the results pertaining to complex-balanced systems, 
for example those in [161 HH [23] , apply to this (deficiency zero and weakly reversible) 
large class of networks. 

It was conjectured at least as early as 1974 that complex-balanced equilibria of 
reaction networks are globally asymptotically stable relative to the interior of their 
positive compatibility classes [35] . This problem was given the name "Global Attrac- 
tor Conjecture" by Craciun et al. [9] , and is considered to be one of the most important 
open problems in the field of chemical reaction network theory [U [S] [Q] [101 [SI] ■ 

Global Attractor Conjecture. A complex-balanced equilibrium contained in the 
interior of a positive compatibility class is a global attractor of the interior of that 
positive class. 

Using the Lyapunov function of Horn and Jackson it is relatively straightforward 
to show that each trajectory of a complex-balanced system remains bounded and 
converges either to the unique equilibrium within the interior of its invariant manifold, 
or to the boundary of the positive orthant, SM^q. Therefore, the Global Attractor 
Conjecture will be proven if it can be shown that any complex-balanced system is 
persistent in the sense of Definition 11.11 below. 

Definition 1.1. For t > denoting time, let zq) be a trajectory to a dy- 
namical system in with initial condition Xq- A trajectory (p{t,Xo) with state space 
R>Q is said to be persistent if 

liminf a;o) > 0, 

for all i £ {l,...,iV}, where (j)i(t,XQ) denotes the ith component of (f>{t,xo). A dy- 
namical system is said to be persistent if each trajectory with positive initial condition 
is persistent. 

We will use the notation (f>{t,xo) for trajectories throughout the paper. We 
see in Definition 11.11 that persistence corresponds to a non-extinction requirement. 
Some authors refer to dynamical systems satisfying the above condition as strongly 
persistent [33j . In their work, persistence only requires the weaker condition that 
limsupj_joo (f>i{t, cco) > for each i e {1, . . . , N}. 

Definition 1.2. Fort > 0, let (j){t,xt^) be a trajectory to a dynamical system in 
with initial condition .tq . The set of w-limit points for this trajectory is the set of 
accumulation points: 

uj{(t){-, xq)) = {x e R"'^ : xq) — > X for some sequence i„ — oo}. 

Note that for bounded trajectories, persistence is equivalent to the condition that 
w(0(-,xo))naM^o = 0. 

It can be shown that a complex-balanced network is necessarily weakly reversible 
[m [131 [20] . Therefore, in light of the above discussion, the Global Attractor Con- 
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jecture is implied by the following, more general, conjecture of Feinberg (see Remark 
6.1.E in [H]): 



Persistence Conjecture. Any weakly reversible reaction network with mass-action 
kinetics and bounded trajectories is persistent. 

Other formulations of the Persistence Conjecture leave out the assumption of bounded 
trajectories, and the above is, therefore, a weaker version of the usual conjecture. In 
fact, it is an open problem as to whether or not weakly reversible networks give rise 
to only bounded trajectories, and we feel it is best to separate these two conjectures. 
Note that the Persistence Conjecture makes no assumption on the choice of rate 
constants. 

Both conjectures remain open. However, in recent years there has been much ac- 
tivity aimed at their resolution. It is known that only certain faces of the boundaries of 
the invariant manifolds can admit w-limit points, those associated with a semi-locking 
set (using the terminology of [1]), which is a subset of the species whose absence is 
forward invariant. (Semi-locking sets were termed siphons in the earlier paper [5] 
in which their concept was formally introduced. However, see Proposition 5.3.1 and 
Remark 6.1.E of [M] for an earlier treatment that anticipated these definitions.) This 
fact has typically focused attention on understanding the behavior of these systems 
near different faces of the boundaries of the invariant manifolds. For example, Ander- 
son [1] and Craciun, Dickenstein, Shiu, and Sturmfels [9] used diflFerent methods to 
independently conclude that vertices of the positive compatibility classes (which are 
polyhedra, see [3]) can not be w-limit points even if they are associated with a semi- 
locking set. In |!3|, it was shown that weak reversibility of the network guarantees that 
facets — faces of one dimension less than the compatibility class itself; that is, a face 
of codimension one — of the positive classes "repel," in a certain sense, trajectories. 
This fact was used to prove the Global Attractor Conjecture when the stoichiomet- 
ric compatibility classes, or invariant manifolds, are two-dimensional. More recently, 
Craciun, Pantea, and Nazarov proved that two-species, weakly reversible systems are 
both persistent and permanent (trajectories eventually enter a fixed, compact subset 
of the strictly positive orthant ]R>o)- They then used this fact to prove that the 
Global Attractor Conjecture holds for three-species systems [10]. Pantea then ex- 
tended these ideas to prove the Global Attractor Conjecture for all systems for which 
the stoichiometric compatibility class has dimension less than or equal to three [29] . 

In [?! the authors studied persistence by introducing the notion of dynamic non- 
emptiability for semi-locking sets, which corresponds to a dominance ordering of the 
monomials near a given face of the compatibility class. These ideas were expanded 
in [25) where the concept of weak dynamic non-emptiahility was introduced, and a 
connection was made to the work on facets in [3] . This work should also be compared 
to the use of strata in both [24J and [51 , where monomial dominance is again considered 
near faces of the invariant manifold. Later, we will see that monomial dominance is 
at the heart of the current paper as well. However, and importantly, the dominance 
is no longer sequestered to individual faces of the invariant manifold and is instead 
considered along sequences of trajectory points in time. Further, the monomials are 
grouped into classes of comparable growth, which allows for a greater understanding 
of the behavior of the system. 

Biological models in which the non-existence of cj-limit points on the boundary 
implies global convergence include the ligand-receptor-antagonist-trap model of G. 
Gnacadja et al. [TSKI^, the enzymatic mechanism of D. Siegel and D. MacLean [50] . 
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and T. McKeithan's T-cell signal transduction model [28 (the mathematical analysis 
appears in the work of E. Sontag [31] and Section 7.1 in the Ph.D. thesis of M. Chavez 
[5]). Recently, Gopalkrishnan showed that any network that violates the Persistence 
Conjecture must be 'catalytic' in a precise sense |19) . 

Complex-balanced systems also play an important role in the study of stochasti- 
cally modeled reaction networks. In [2^, Anderson, Craciun, and Kurtz showed that 
a stochastically modeled system admits a product form stationary distribution if the 
associated deterministically modeled system admits a complex-balanced equilibrium 
(see also [271 in which the same results were derived independently by David Luben- 
sky). Also, in the study of stochastically modeled reaction systems with multiple 
scales, which is the norm as opposed to the exception in the stochastic setting, it is 
often desirable to perform a reduction to the system using asymptotic analysis (usu- 
ally averaging and law of large number techniques), see, for example, [7] or [26| . If one 
component of the system can be shown to behave deterministically in the asymptotic 
limit, then knowing this component converges to a steady state (i.e. knowing that the 
conclusions of the current paper hold) may allow for the proof of results pertaining 
to the dynamics of the other components. 

1.2. Results in this paper. In this paper, we will provide multiple results 
pertaining to deterministically modeled, weakly reversible systems. Most results will 
pertain to the dynamics of individual trajectories. These results will allow us to 
conclude that the Global Attractor Conjecture holds in the case when the underlying 
reaction network consists of one linkage class, or connected component. It is worth 
noting that we will not provide a proof of the Persistence Conjecture in the one linkage 
class case. As will become apparent, the technical difference between the conjectures 
will be captured by a condition on where the w-limit points of a trajectory can reside, 
see Theorem 14. 101 

To prove our results, we will introduce a method for partitioning the relevant 
monomials of the dynamical system along sequences of trajectory points into classes 
with comparable growths. This method will allow us to conclude that if a trajectory 
converges to the boundary, then a whole family of Lyapunov functions decrease along 
the trajectory. We will then be able to overcome the fact that the usual Lyapunov 
functions of chemical reaction network theory are bounded on the boundary of the 
positive orthant, which has been the technical sticking point to a proof of the Global 
Attractor Conjecture in the past. The methods developed should prove useful in 
future contexts, both deterministic and stochastic, as well as the current one. Also, 
it will be natural to focus our attention on systems with a generalized mass-action 
kinetics in which the rate constants are allowed to be functions of time. This context 
is useful because the projection of a trajectory of a reaction system onto some relevant 
subset of the species can itself be viewed as a trajectory of a reaction system with 
generalized mass-action kinetics. 

The outline of the paper is as follows. In Section |2l we will provide the requisite 
definitions and terminology from chemical reaction network theory. In Section [31 we 
will discuss projected dynamics, and introduce and develop the basic properties of 
reduced reaction networks and generalized mass-action systems. In Section HJ we will 
give our main results together with their proofs. In Section jSj we present an example 
to demonstrate our results. 

2. Preliminary concepts and definitions. Most of the following definitions 
are standard in chemical reaction network theory. The interested reader should see 
[T3] or [50] for a more detailed introduction. 
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Reaction networks. An example of a chemical reaction is 25*1 + 5*2 ^3, where we 
interpret the above as saying two molecules of type 5*1 combine with a molecule of type 
S'2 to produce a molecule of type 5*3. For now, assume that there are no other reactions 
under consideration. The Si are called chemical species and the linear combinations 
of the species found at either end of the reaction arrow, namely 25*1 + S'2 and ^3, are 
called chemical complexes. Assigning the source (or reactant) complex 2Si + S'2 to 
the vector y = (2, 1, 0) and the product complex S3 to the vector y' = (0, 0, 1), we can 
formally write the reaction as y — t- y'. 

In the general setting we denote the number of species by N, and for i G 
{1, . . . ,N} we denote the ith species as Si. We then consider a finite set of reac- 
tions with the fcth denoted by yk y'f^, where yk,yk £ ^>o ^■'"S (non-equal) vectors 
whose components give the coefficients of the source and product complexes, respec- 
tively. Using a slight abuse of notation, we will also refer to the vectors yk and 
as the complexes. Note that if yj, = or = for some k, then the fcth reaction 
represents an input or output, respectively, to the system. Note also that any complex 
may appear as both a source complex and a product complex in the system. We will 
usually, though not always (for example, see condition 3 in Definition 12.11 below) use 
the prime ' to denote the product complex of a given reaction. 

As an example, suppose that the entire system consists of the two species Si and 
S2 and the two reactions 

Si ^ S'2 and S2 ^ Si, (2.1) 

where Si S2 is arbitrarily labeled as "reaction 1." Then N = 2 and 

yi = (l,0), 2/'i = (0,l) and 2/2 -(0,1), 2/2 = (l,0). 

Thus, the vector (1,0), or equivalently the complex Si, is both j/i, the source of the 
first reaction, and 2/2: the product of the second. 

For ease of notation, when there is no need for enumeration we will typically drop 
the subscript k from the notation for the complexes and reactions. 

Definition 2.1. Let S = {S^jfL-^^, C = {y} with y e Z^^, and U = {y ^ 
y'} denote finite sets of species, complexes, and reactions, respectively. The triple 
{S,C,TZ} is called a chemical reaction network so long as the following three natural 
requirements are met: 

1. For each Si G S, there exists at least one complex y €z C for which yi > 1. 

2. There is no trivial reaction y ~> y € TZ for some complex y & C. 

3. For any y £ C, there must exist ay' £ C for which y ^ y' £ TZ or y' ^ y G TZ. 
If the triple {S,C,TZ} satisfies all of the above requirements except 1., above, then we 
say {S,C,TZ} is a chemical reaction network with inactive species. 

Notation: We will use each of the following choices of notation to denote a 
complex from C: y, y\ yk, y'^, Vi, Vj, ye, and even z^-. However, there will be other 
times in which we wish to denote the ith component of a complex. If the complex 
in question has been denoted by yk, then we will write ykj. However, if the complex 
has been denoted by y, then we would write its ith component as yi, which, through 
context, should not cause confusion with a choice of complex yi. See, for example, 
condition 1 in Definition 12.11 above. 

Definition 2.2. To each reaction network {S,C,Tl} we assign a unique directed 
graph called a reaction diagram constructed in the following manner. The nodes of the 
graph are the complexes, C. A directed edge {y,y') exists if and only if y y' £ TZ. 
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Each connected component of the resulting graph is termed a linkage class of the graph. 



For example, the system described in and around (j2.ip has reaction diagram 
Si ^ S2, which consists of a single linkage class. 

Definition 2.3. Let {S,C,TZ} denote a chemical reaction network. Denote the 
complexes of the ith linkage class by Li C C. We say a T <Z C consists of a imion of 
linkage classes if T — Un^iLi for some nonempty index set I. 

Definition 2.4. The chemical reaction network {S,C,TZ} is said to be weakly 
reversible if each linkage class of the corresponding reaction diagram is strongly con- 
nected. A network is said to be reversible if y' ^ y G TZ whenever y ^ y' G TZ. 

It is easy to see that a chemical reaction network is weakly reversible if and only 
if for each reaction y ^ y' ^ TZ, there exists a sequence of complexes, . . . , j/^ G C, 
such that y' — > yi G 7^, t/i — > 2/2 G 7?., • • ■ , Ur-i TJr ^TZ, and j/,. y eTZ. 

Dynamics. A chemical reaction network gives rise to a dynamical system by way 
of a rate function for each reaction. That is, for each y^ y'^. € TZ, or simply 
k £ {1, . . . , |7^|}, we suppose the existence of a function Rk = Ry^^y'^ that deter- 
mines the rate of that reaction. The functions Rk are typically referred to as the 
kinetics of the system and will be denoted by /C, or lC(t) in the non-autonomous case. 
The dynamics of the system is then given by the following coupled set of (typically 
nonlinear) ordinary differential equations 



where k enumerates over the reactions and x{t) G M>o ^ vector whose ith component 
represents the concentration of species Si at time t>Q. 

Definition 2.5. A chemical reaction network {S,C,TZ} together with a choice 
of kinetics IC is called a chemical reaction system and is denoted via the quadruple 
{S,C,TZ,]C}. In the non- autonomous case where the Rk can depend explicitly on t, we 
will write {S, C, TZ, JC(t)}. We say that a chemical reaction system is weakly reversible 
if its underlying network is. 

Integrating (|2.2p with respect to time yields 



Therefore, x(t) — x(0) remains within S — span{y^ — yk}ke{i....,R} for time. 

Definition 2.6. The stoichiometric subspace of a network is the linear space 
S — span{y^ — yk}k£{i,...,\n\}- The vectors y'f. — yk are called the reaction vectors. 

Under mild conditions on the rate functions of a system, a trajectory x{t) with 
strictly positive initial condition a;(0) G M^^q remains in the strictly positive orthant 
M^g for all time (see, for example. Lemma 2.1 of [31]). Thus, the trajectory remains 
in the relatively open set {x{0) + S) n R^q, where x(0) + S -.^ {z € \ z = 
x{0) + V, for some v G S}, for all time. In other words, this set is forward-invariant 
with respect to the dynamics. It is also easy to show that under the same mild 
conditions on Rk, {x{0) -f S) nM>g is forward invariant with respect to the dynamics. 
The sets (a;(0)-f S')nR^Q will be referred to as the positive stoichiometric compatibility 
classes, or simply as the positive classes. 




(2.2) 



k 
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The most common kinetics is that of mass- action kinetics. A chemical reaction 
system is said to have mass-action kinetics if aU rate functions R/. ~ Ry^^y'^. take the 
multiphcative form 

i?fe(:E) = (2.3) 

where is a positive reaction rate constant and yk is the source complex for the 
reaction. For u G K>q and v G M^, we define 

where we have adopted the convention that O'' = 1 , and the above is undefined if = 
when Vi < 0. Mass action kinetics can then be written succinctly as Rk{x) — k^x^'' . 
Combining (|2.2p and (|2.3p gives the following system of differential equations, which 
is the main object of study in this paper, 

i(i) = ^«;,a;(t)«'=(y;: -2/fe). (2.4) 

k 

While it is the properties of solutions to the system (|2.4p that are of most interest 
to us, it will be natural for us to consider systems with a generalized form of mass- 
action kinetics. The following definition is similar to Definition 2.6 in |10| for "k- 
variable mass-action systems." See also [6j for a recent treatment of chemical reaction 
systems with non-autonomous dynamics. 

Definition 2.7. We say that the non- autonomous system {S,C,TZ,JC{t)} has 
bounded mass-action kinetics if there exists anrj > such that for each k {I, . . . ,\TZ\} 

Rk{x,t) ^ Kk(t)xy\ 

where rj < Kfe(t) < l/rj for all t > and k e {1, . . . , \TZ\}. 

2.1. Complex balanced equilibria and the deficiency of a network. The 

Global Attractor Conjecture, which was stated in Section 11.11 is concerned with the 
asymptotic stability of complex-balanced equilibria for mass-action systems. For each 
complex y e C we will write {k \ = y} and {k \ = y} for the subsets of reactions 
Vk ^ y'k ^ for which y is the source and product complex, respectively. 

Definition 2.8. We say c is an equilibrium of the dynamical system x{t) = 
f{x{t)), if f{c) — 0. An equilibrium c G M>q of (|2.4p is said to be complex-balanced 
if the following equality holds for each complex y E C: 

^ KkC^" = ^ KfcC^'". 

{fc I yk=y} {k I y'k=y} 

Note that on the right hand side of the above equality, yk represents the source 
complex for a given reaction for which y is the product complex, whereas on the left 
hand side each source complex is identically equal to y. Thus, c is a complex-balanced 
equilibrium if for all complexes y G C, at concentration c the sum of rates for reactions 
for which y is the source is equal to the sum of rates for reactions for which y is the 
product. That is, crudely, the total flux into complex y is equal to the total flux out 
of complex y. A complex-balanced system is a mass-action system {S,C,TZ,JC} that 
admits a complex-balanced equilibrium with strictly positive components. 
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In [5], complex-balanced systems are called "toric dynamical systems" in order 
to highlight their inherent algebraic structure. There are two important special cases 
of complex-balanced systems: the detailed-balanced systems and the zero deficiency 
systems. 

Definition 2.9. An equilibrium c G R>q of a reversible system with dynamics 
given by mass-action kinetics (|2.4|) is said to be detailed-balanced if for any pair of 
reversible reactions yk ^ y'^ with forward reaction rate Kk and backward rate n'^ the 
following equality holds: KkC^'' — k'^.c^''. 

A detailed-balanced system is a reversible system with dynamics given by mass- 
action kinetics (I2.4p that admits a strictly positive detailed-balanced equilibrium. It 
is clear that detailed-balancing implies complex-balancing. 

Definition 2.10. For a chemical reaction network {iS,C,7?,}, let n denote the 
number of complexes, I the number of linkage classes, and s the dimension of the 
stoichiometric subspace, S. The deficiency of the reaction network is the integer 
n — I — s. 

The deficiency of a reaction network is non-negative because it can be interpreted 
as either the dimension of a certain linear subspace [13] or the codimension of a 
certain ideal [9]. Note that the deficiency depends only on the reaction network and 
not the choice of kinetics. The Deficiency Zero Theorem tells us that any weakly 
reversible reaction system (12. 4p whose deficiency is zero is complex-balanced, and 
that this fact is independent of the choice of rate constants Kk [H]. On the other 
hand, a reaction diagram with a deficiency that is positive may give rise to a system 
that is both complex- and detailed-balanced, complex- but not detailed-balanced, or 
neither, depending on the values of the rate constants Kk [9] fT2 l [TS l [2T|. 

Complex-balanced systems have the property that there is a unique, complex- 
balanced equilibrium within the interior of each positive stoichiometric compatibility 
class [m [22j [23] . Thus, proving that each such equilibrium is globally asymptotically 
stable relative to its positive class, i.e. showing the conclusion of the Global Attrac- 
tor Conjecture holds, would completely characterize the long-time behavior of these 
systems. 

3. Projected dynamical systems and reduced reaction netv^rorks. The 

two related concepts presented below, projected dynamics and reduced reaction net- 
works, will be used in the proofs of the main theorems of Section [4l There, we will 
consider a bounded trajectory of a system and incorporate the dynamics of those 
species which do not approach the boundary of the positive class into the rate con- 
stants, thereby yielding a system with bounded mass- action kinetics. This will be the 
projection of the dynamics onto those species that do approach the boundary. The 
resulting reaction network will be the reduced network. 

3.1. Projected dynamics. As discussed in Section [21 our interest is in the 
qualitative dynamics of an N dimensional, autonomous systems of differential equa- 
tions with parameters k = . . . , K|7?,|). That is, we are considering systems of the 
general form 

il{t) = fi{K,Xi{t), . . .,XNit)) 

: (3.1) 

XNit) = fN{K,Xi{t), . . .,XN{t)). 

To study these systems, it will be natural to later consider an associated non- 
autonomous set of differential equations constructed by projecting (13.11) onto a subset 
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of the dependent variables. Specifically, let J7 C {1, . . . , N} be nonempty with \U\ = 
M < N. Without loss of generality, assume for now that U = {1, . . . , M}. We now 
consider the dynamical system in R*^ with state vector x\ij, whose ith component, 
denoted x\u,i, is for all time equal to Xi, defined in (|3.ip . We see that the vector 
valued function x\u satisfies 

x\u,i{t) = fl{K,x\u,l{t), ■ ■ . ,x\u,Mit),XM+l{t), ■ ■ . ,XAr(i)) 
= fl{K{t),x\u,l{t), ■ . .,x\u,M{t)) 

: (3.2) 

a;|(7,M(i) = fMiK,x\u,l{t), . . . ,x\u,Mit),XM+l{t), ■ ■ .,XN{t)) 
= fM(Mt)iX\uA{t), ■ ■ ■,x\u.M{t)), 

where k(t) = {k, xm+i (i), • ■ • , a;jv(t)) and fi, i G {1, . . . , M}, are defined via the above 
equalities. Thus, the system x\u.i satisfies the same differential equations as does the 
system for Xi, except any dependence upon the variables Xj, for j > M + 1, has been 
incorporated into the dynamics as known, time-dependent functions. That is, they 
are now viewed as inputs, or perhaps forcing, to the system. We will call the system 
p.2|) the projected dynamics of p.f p with respect to U. 
For example, consider the system 

xi ~ ~nixix\ — K2X1X3 + K5X2 

±2 = K^x^ - 2kiXixI - ^5X2 , (3.3) 
is = K4 + Kixix\ — Ka^a 

where k = {ki, . . . , K5) € R>o- Then for U = {1, 3} the projected dynamics of p.3p 
with respect to U is 

Xl = -HiiC,l(t)xi - K2X1X3 + K5C2(i) 
is = '«4 + KlCl{t)xi - K3X3 

where (i{t) — X2{t)'^, C,2{t) ~ X2{t), and X2{t) is still defined via the system p.3p . 
The goal would now be to translate any control we can get over X2{t), and hence C(^)) 
into qualitative information about the behavior of xi and x^. Later, we will simply 
incorporate the function C of p.4[) into the variables Kfc and view each Kk{t) as a 
function of time. 

3.2. Reduced reaction networks. We begin with more notation. Let v £ 
for some iV > 1, and let U C {1, . . . , N} be nonempty. We write U[j] for the jth. 
component of U. We then write v\u to denote the vector of size \U\ with 

V\u,j = iv\u)j =VU[3] 

for j G {I, . ■ . ,\U\}. Thus, v\ij simply denotes the projection of v onto the components 
enumerated by U. For example, if iV = 8 and U — {2,4,7}, then for any v G R^, 

v\u = {V2,V4,V7). 

Definition 3.f. Consider a reaction network {S,C,TZ} with S = {Si, . . . , Sn} 
and let U C {1, . . . , N} be nonempty. The reduced reaction network of {S,C,TZ} 
associated with U is the reaction network {Su,Cu,TZu} constructed in the following 
manner: 
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1. SetSu = {S'j e 5 : i^U). 

2. Set Cjj = {y\ij : y £ C}. We say the complex y reduces to the complex y\ij- 

3. SetTZu = {y\u ^ y\'u ■ y^y'^T^i and for which y\u ^ y'\u}- 

4- If a resulting linkage class consists of a single complex, we delete that complex 
from Cu . 

Example 3.2. Consider the reaction network with species S = {Si, . . . , S^} and 
reaction diagram 

2Si + S2^S3 + 2Si 

S's + 5*3 <^ 5*1 + S'4 <^ S's + <S'2 

K7 

5*5 + 25*2 <=^ 52 + 53, 

Kg 

where we have ordered the reactions and incorporated the rate constants into the re- 
action diagram. Let U — {1,4,5}. Then, Su = {Si, S4, S5}, Cjj = {Sz,Si + 5*4,0}, 
and the resulting diagram of the reduced reaction network is 

0^55^5i + 54. (3.6) 

Here, both the complex 25i + 52 and S3 + 25i reduce to 25i. However, by rule 3 in 
Definition \3.1\ we do not include 25i — >■ 25i in TZu, and by rule 4 we delete 25i from 
Cjj- Note also that the original network has three linkage classes whereas the reduced 
network has only one. 

Note that because of rule 4 in Definition 13. 1[ it is possible to have Si € Sjj, 
even though Si does not appear in any complex in Cu- For example, if 1 € C/, but 
2,3 ^ U, and the only reactions in which 5i participates are 5i + 52 ^ 5i + 53, 
then 5i G Su, even though Si does not appear in any complex in Cu- In this case, 
the reduced reaction network has inactive species, see Definition 12.11 Note, however, 
that this situation only arises if the concentration of 5i was time independent in the 
original system. Thus, the original system could have been reduced by incorporating 
Si into the rate constants. Such an incorporation can have the effect of lowering the 
deficiency of the network without changing the dynamics (see [11] to see a treatment 
of how different network structures may give rise to the same dynamical system). 

The following lemmas give some insight into how the structure of {Su,Cu,Ti-u} 
depends upon the structure of {S,C,TZ}. Both will be used in Section|4]in the proof 
of our main results. 

Lemma 3.3. Let {S,C,TZ} be a reaction network with S — {5i, . . . , 5jv}. Let 
U C {1,...,A^} be nonempty. Then, the reduced reaction network {Su,Cu,Ti-u} has 
less than or equal to the number of linkage classes as {S,C,TZ}. 

Proof. Condition 3 of Definition 13 . 1 1 shows that if yi , 2/2 G C are in the same linkage 
class, then yi\u and 2/2 Ic/ are also. Thus the result is shown simply by counting the 
number of unique linkage classes of {S,C,TZ} and {Su,Cu,TZu} by enumerating over 
the complexes C and Cu, respectively. □ 

Lemma 3.4. Suppose that {S, C, TZ}, with S = {Si, . . . , Sn}, is weakly reversible 
and that U C {1, . . . , N} is nonempty. Then {Su,Cu, Ti-u} is weakly reversible. 

Proof. Suppose y\u y\'u G 72.. By construction there are complexes y,y' € C, 
with y ^ y' & TZ, that reduce to y\u,y\'u- By the weak reversibility of {5, C, TZ}, there 
is a sequence of directed reactions, yk — > G 72, beginning with y' and ending with y. 
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(3.5) 



If for each — ^ of this sequence we have yk\u y'k\u i then, by construction, there 
is a sequence of directed reactions in TZjj beginning with y\'jj and ending with y\ij. If 
yk\u = y'k\u for one of the k, such a sequence still exists except now yk\u y'kj^i\u- 
Similar reasoning holds when yk\u = y'k\u for more than one k. □ 

We need to provide the reduced network {Sjj ^Cjj ^Ti-u} with a natural kinetics. 
The kinetics, lC{t), is given via the projection of the dynamics, described in Section 
13.11 of {5,C,7^, /C} onto the elements of U. The variables Hk{t) are now functions of 
time. Note that the dynamics of the resulting projected system depends upon the 
dynamics of the original system. 

Example 3.5. Again consider the reaction system with species S — {5*1, . . . , 6*5} 
and reaction diagram (|3.5|) . For U = {1,4,5}, the reduced network was p.6|) . Incor- 
porating the projected dynamics yields 

^Ki(t) K3(t) 
<^ ^5 S'l + ^4, 

K2(t) ki(t) 

where ki{t) = KiiX2{t)xz{t), k2{t) = K^X2{tY, k^{t) = K,?,x^{t) + KeX2{t), ki{t) = 
(K4 + K5), and the indexing k of the kk{t) is over reactions in the reduced network and 
not the original network. 

It is important to note that the variables kk (t) for the reduced system, which take 
the place of the rate constants, are always non-negative as they consist of positive 
linear combinations of non-negative monomials of the variables Xj for which j ^ 
U, see p.7p below. Also, while the reduced system is a non-autonomous system 
with generalized mass-action kinetics, the functions kk{t) are not necessarily bounded 
either above or below. Finally, note that the functions kk{t) depend explicitly on the 
original trajectory of the system and, in particular, will be different functions of time 
for different initial conditions of the system {S,C,TZ, K.}. 

It is useful to have a more formal representation of the projected dynamics. Thus, 
let {S^C,TZ,JC} be a reaction network with mass-action kinetics, K. = {nk}- Let 
[/ C {1, . . . , be nonempty. The reduced mass-action system of {5, C, TZ, K.} with 
respect to U is the non- autonomous mass-action system {Su,Cuj'R-u,l^u(t)}, with 
^u{i) = where the indexing k of the K,k{i) is over reactions in the reduced 

network and not the original network, and where for yk\u — ^ y'k\u 6 T^u 

^kit)^ K,ixit)\u^r^"\ (3.7) 

{zi^z'-eTZ : yk\u=Zi\u and j,^|u=zj|[/} 

where x{t) is the solution to equation (|2.4p for the system {S,C,TZ, K.}. 

We reiterate the fact that based upon the above definitions, the differential equa- 
tions governing the dynamics of Xi for i G U for the reduced system are exactly the 
same as the differential equations for Xi for i G [/ of the original mass-action system. 
In Section 21 we will project the dynamics of a bounded trajectory onto the subset 
of the species that go to the boundary, thereby producing a trajectory that satisfies 
a system with bounded mass-action kinetics. Note that both the resulting reduced 
network and the bounds on the kinetics will depend upon the initial condition of the 
original system. 

4. Main results. We begin in Section |4T] bv introducing the concept of parti- 
tioning a set of vectors along a sequence, which will be one of our main analytical 
tools and will allow us to group the relevant monomials of the dynamical system 
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along sequences of trajectory points into classes with comparable growths. In Sec- 
tion 14.21 our main results will be stated and proved. The main results will focus on 
non- autonomous systems because we will later project our system of interest, which is 
autonomous, to the non- autonomous system consisting of those species that approach 
the boundary along a subsequence of times. 

4.1. Partitioning vectors along a sequence. We begin by recalling that for 
any vectors u, v such that u G R>q and v € we define u" = u^^ ■ • ■ , where 
we use the convention O" = 1. All sequences and subsequences will be indexed by 
non-negative integers. 

Definition 4.1. LetC denote a finite set of vectors mM^. Let Xn & K>o denote 
a sequence of points in the strictly positive orthant. We say that C is partitioned along 
the sequence {a;„} if there exists Ti C C, i — 1^ . . . , P, termed tiers, such that Ti ^ 0, 
TiClTj = if i ^ j , and UiTi — C (that is, the tiers constitute a partition of C), and 
a constant C > 1, such that 
(i) */ J/j J Vk G Ti for some i G {1, . . . , P}, then for all n 

}_T,Vj < ^Vk < n^Vj 
which is equivalent to either of the conditions 

^<^<C or -<xl--y^<C. 
(ii) if yk G a,nd yj G T^+m for some m G {1, . . . , P ~ i} , then 

n 

— ^- — > oo, as n —i> oo. 

Xn 

Therefore, we have a natural ordering of the tiers: Ti )^ T2 >- T3 >-■••>- Tp, and 
we say Ti is the "highest" tier, whereas Tp is the "lowest" tier. 

Throughout the paper. Definition 14.11 will be used in the context of the sequence 
{xn} being trajectory points, and the set of vectors C being the complex vectors. 

The following lemma, which is critical for our purposes, states that given a set 
of vectors and a sequence of points in R^q, there always exists a subsequence along 
which the vectors are partitioned. 

Lemma 4.2. Let C denote a finite set of vectors in M^. Let Xn he a sequence of 
points in M>q. Then, there exists a subsequence of {x„} along which C is partitioned. 

Proof. Denote the elements of C as 1 < i < r, where |C| = r. Note that there 
are r! < 00 ways to order the elements of C. Therefore, there exists a reordering of 
the vectors in C such that the set of indices n for which 

xl->xl^>---> xl (4.1) 

holds is infinite; letting denote those indices, we obtain the corresponding subse- 
quence {xji^}. Thus, we have instituted an ordering, y\ >- y2 ■ ■ ■ ^ Vr along this 
subsequence, and y\ can be viewed as maximal. The goal now is to simply get more 
information about this ordering (along further subsequences) and ask which vectors 
stay "close" to each other, and which diverge. This will give us the natural dividing 
lines for our tiers. 



12 



For i £ {l,...,r - 1}, define -0* : ^>o K by ipi{x) = x^^jx^^^K By the 
inequalities (14.11) . for each i G {1, . . . , r — 1} and for fc > 1, we have ^iixn^) > 1- We 
will construct the tiers. We begin by setting T\ — {yi}. Next, we ask if 

lim inf -01 (a;„j, ) < oo. (4-2) 

If (14. 2p holds, we set T\ — {2/1,2/2} and redefine our sequence {a^n^,} as an appropriate 
subsequence so that linifc-^oo ''Pii^rik) exists, and is finite. Next, we ask if 

lim inf -02 (a;nfc ) < 00 

along this new subsequence. If so, we set Ti = {2/1,2/2, 2/3} and redefine our sequence 
appropriately so that linifc_^oo "02 {xn^ ) exists and is finite. This process will either 
terminate with Ti = C or when liniinffe_j.oo i'b-iixn,,) = 00 for some b G {2, . . . , r}. If 
such a b exists, we have Ti = {2/1, . . . ,2/6-1}, and then we begin building the second 
tier by setting T2 = {yt}- Now fill T2 in the same manner that we did Ti by looking at 
the values of liminffc_j.oo ^i{xnk) foi' the appropriate i's. Repeat this process, always 
redefining the sequence as the subsequence guaranteed to exist at each step, until we 
have a sequence of tiers Ti, T2, . . . , Tp. 

It remains to find the appropriate C > 1 for Definition 14.11 For each pair 
2/j,2/j+i G Ti, we define C^j- == limfc_^oo a^ni/a^ni^^ which exists and is finite by con- 
struction. For each such i, j, let Cij — Cij + 1 and note that there is a if > 1 so 
that 

1 < < r 

for all k > K, and all relevant i, j (that is, this bound is uniform in i and j). Let 

C = max TT Ci j . 

i6{l,...,P} 

{j ■■ y],V]+i&Ti} 

Then, restricting ourselves to the subsequence with k > K , C is now partitioned along 
the resulting subsequence with tiers {Ti}, i — 1, . . . , P, and constant C > 1. □ 

The following lemma states that for any set of vectors in M™, either their span 
includes a non-zero vector in the non-positive orthant IR™q, or there is vector normal 
to their span that intersects the strictly positive orthant. 

Lemma 4.3 (Stiemke's Theorem, [32]). For i ^ 1, . . . ,n, let m e M™. Either 
there exists an a £ M" such that 



. , TO 



and such that at least one of the inequalities is strict, or there is a w E R"q such that 
w ■ Ui — for each i G {1, . . . , n}. 

Definition 4.4. Let w eR'^ . The set {i e {1, . . . ,N} : Wi ^ 0} is called the 
support of w. 

Definition 4.5. Let C denote a finite set of vectors in R^. Let {Ti} denote a 
partition of C. Let U C {1, . . . ,N} be nonempty. We say that the vector w G M>q is 
a non-negative conservation relation that respects the pair (U, {Ti}) if the following 
two conditions hold: 
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1. Wi > if and only if i ^ U . That is, the support of w is U. 

2. Whenever yj, yg £ Ti for some i, we have that w ■ {yj — yi) = 0. 

Note that if each Ti consists of a single element, then any vector w whose support 
is U satisfies the requirements of the definition as the second condition holds auto- 
matically. Also note that if C = Ti , then the type of conservation relation described 
above is the usual concept in chemical reaction network theory. 

If in the following theorem, C is taken to be the set of complexes of a reaction 
network, then the theorem guarantees that there must be a conservation relation 
among certain subsets of the reactions if a trajectory converges to the boundary of 
the positive orthant. 



Theorem 4.6. Let C denote a finite set of vectors in M.^ . Let a;„ € M^q denote 



a sequence of points with Xn ^ z € dM.^^, as n oo. Let U = U{z) — {i G 
{1, . . . ,N} : Zi = 0}. Finally, suppose that C is partitioned along {a;„} with tiers Ti, 
for i ~ 1,...,P, and constant C > 0. Then, there is a non-negative conservation 
relation w G M>g that respects the pair {U, {Ti}). 

Proof. We suppose, in order to find a contradiction, that there is no non-negative 
conservation relation that respects the pair ([/, {T;}). Define the sets Wi C M^, for 
i = 1, . . .P, and C via 

p 

= {yj - yk I yj,yk G TJ, W = |J W,, 

1=1 

and denote the elements of W by {ukj- Note that if Ti consists of a single element, 
then Wi consists solely of the vector 0. Let m — \U\ > be the number of elements 
in U and let Wi\u C K™ and W\u C M" be the restrictions of Wi and W to the 
components associated with the index set U, as discussed in Section Denote the 
elements of Ty|(7 by {vk}- Thus, collecting terminology, Uk € K^, whereas Vk G K"*, 
and for each Uk G W, there is a corresponding £ W\{j for which Uk\u = Vk, however 
the mapping -lu need not be injective. 

The set Wlu must contain at least one nonzero vector because otherwise any 
non-negative vector with support U would be a non-negative conservation relation 
that respects the pair {U, {Ti}), but we have assumed that no such relation exists. 

Because we have assumed there is no conservation relation that respects the pair 
{U,{Ti}), we may conclude by Lemma [43l that span(Vl^|(7) must intersect IR™g in a 
non-trivial manner. That is, there exist c^^ G M such that 



E 



CkVk < 0, (4.3) 



for each j € {1,...,to}, and such that the inequality is strict for at least one j. 

For Vk G W\u, let denote the number of vectors of W that reduce to it (recall 
that the operation need not be injective). Define the function M : R>q ^ R by 



M{x)'= n 



where Ck and are chosen for Uk G if Uk\u — Vk G W\u. Note that, by 
construction and by the definition of partitioning along a sequence, if Uk G W , then 
there are yj,yi G Ti for some i, such that Uk — ye — Uj and 

1 j-Vi 

_ < ™"fc _ ZJL. < C 
/-< — " ~ — ^' 
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for all n > 1. Therefore, the sequence M(a;„) is uniformly, in n, bounded both from 
above and below. Noting that each Xn has strictly positive components, we may take 
logarithms and find 

ln(M(x„)) = ( —Uk]-lnxn, 

where for a vector u G M^o define 

ln(u) (In(wi), • • • , In('UAr)). 
Expanding along elements of U and yields, 

ln(M(a;„)) = ( ^ c^Wfe j ■ ln(x„|(7) + f ^ — u^lc/- j • ln(x„|c/0- (4.4) 

By construction, Xn.i is bounded from both above and below for £ d W^. Thus, the 
second term in (j4.4p is bounded from above and below. By the inequality (14. 3p , where 
at least one term is strictly negative, and the fact that Xn.j for each j d U along 
this sequence, we may conclude that the first term, and hence ln(M(a;„)) itself, is 
unbounded towards +oo, as n oo. This is a contradiction with the previously 
found fact that the sequence M{x„) is uniformly bounded above and below, and the 
result is shown. □ 

4.2. Persistence and the Global Attractor Conjecture in the single link- 
age class case. For any x S R>^oi define Vx '■ -K>o ^ ^>q by 

N 

V^{X) = i^^iHx^) - H^^) - 1) + X,] . (4.5) 

i=l 

For X e dR^Q, define V^ix) via the continuous extension of (|4.5p . This is the standard 
Lyapunov function of chemical reaction network theory [T31 [10] and has commonly 
been used in the study of complex-balanced systems where x is typically taken to be 
a complex-balanced equilibrium. However, we emphasize that in the current setting 
X need not be a complex balanced equilibrium. Note that Wx{x) —Inx — Inx. It is 
relatively straightforward to show that for any x G K>o, Vx is convex with a global 
minimum of zero at a; [13| . 

The outline of the technical arguments to come is as follows. In Lemma 14.71 we 
will show that for a trajectory of a non-autonomous, weakly reversible system with 
bounded kinetics to approach the boundary, at least one of two conditions, CI or C2, 
must hold. In Lemma l4.8l we will show that condition C2 can not hold for trajectories 
of systems of interest in this paper. Condition CI of Lemma 14771 will then essentially 
tell us that a whole family of Lyapunov functions decrease along the trajectory. This 
is a substantially stronger condition than having a single Lyapunov function decrease 
along a trajectory, and will eventually allow us, with the help of Lemma 14.91 to 
overcome the fact that Vx{x) does not diverge to -l-oo as x ^R^g, which has been 
the technical sticking point to a proof of the Global Attractor Conjecture in the past. 

Lemma 4.7. Let {S,C,TZ,]C{t)}, with S = {Si, . . . , Sn} , be a weakly reversible, 
non- autonomous mass-action system with bounded kinetics. For t > 0, let x{t) = 
(j){t,xo) be the solution to the system with initial condition xq. Suppose xq G M>o 
such that (l){t,xo) remains bounded and dist((/>(t, xq), 9K>q) — )■ 0, as t oo. Then at 
least one of the following two conditions hold for this trajectory: 
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CI: For any x e IR>0' there exists a T — > such that t > T implies 



dt 



V^{x{t)) - ^k{t)x{t)y- (y^ - yu) ■ {Hx{t)) - \n{x)) < 0. 



OO, such that Xn = (j){tn,XQ) G 



N 
>0 



C2: There exists a sequence of times, tn 

converges to a point z G uj{(j){-, xo)) H 9IR>q, and 
(i) C is partitioned along Xn with tiers {Ti}fLi, and constant C, and 
(ii) Ti consists of a union of linkage classes. 
Proof. We suppose condition CI does not hold, and will conclude that condition 
C2 must then hold. Because condition CI does not hold, there is an a; G M>q and a 
sequence f„ — >■ oo such that 



E 

k 



Kk{tn)xl''{y'k - yk) ■ (ln(a;„) - ln(a;)) > 0, 



(4.6) 



where Xn = 4>{tn, xq). We now fix x. 

Combining dist(0(t„, xq), 9R>o) ^ 0? co, with the boundedness of the 

trajectory allows us to conclude that there exists a convergent subsequence of {x„}, 
which we take to be the sequence itself, with limit point z G a;((/)(-, xq)) n dM.^Q. 
Note that by construction the inequality (14.61) holds for all Xn of the subsequence. 
Applying Lemma 14.21 we partition the complexes along an appropriate subsequence 
of the sequence, which we will again denote {xn}, with tiers T^, i = and 
constant C > 1. 

In the following, for tier i G {1, . . . , P}, we denote by 

• {i — i} all reactions with both source and product complex in T^, 

• {i ^ i + m} all reactions with source complex in Ti and product complex in 
T:,+™ for m G {1, . . . , P - i}, 

• {i ^ i — m} all reactions with source complex in Ti and product complex in 
Ti_„i for TO G {1, . . . , i - 1}. 

Defining u/v = (ui/vi, . . . ,un/v^) for w, u G K>oj may re- write the left hand side 
of the inequality (|4.6p 



i=l {i— yi} 

+ Kkitn)a 
m=l {i-^i+rn} 



i-1 

m=l {i^i-m} 



In 



3fc 



(4.7) 



■ Cfc 



(4.8) 



(4.9) 



where for the kth reaction Ck = ln(a;^''/a;^'») — —{y'l. — yk) - Inx. Note that supj, \ck\ < 
OO because x is fixed. Note also that, by construction, for large enough n any compo- 
nent in the enumeration ()4.8p is negative, and, in fact, \n(xn'' / xl^'') — >■ — oo as n ^ oo, 
for these terms. We will now show that the total summation above (that is, the left 
hand side of (|4.7I) and, hence, (|4.6p ) must also, for large enough n, be strictly negative 
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unless condition C2 holds. This will then conclude the proof as it shows "not CI =^ 
C2." 

Suppose condition C2 does not hold. Then, for the specific partition we have 
along {xn}, it must be that Ti does not consist of a union of linkage classes. Thus, 
by the weak reversibility of the system there must be at least one reaction, yki -> y'/.^ 
say, such that y^i G Ti and y'/^i & Tj for j > 2. That is, there is a reaction being 
enumerated in (14.81) with i = 1 and m > 1. As noted above, we have by construction 

that \n{xn''^ /xfil''-) — — oo, as n — > oo. Further, there is a C > 1 such that for any 

yj G C, xl^''^ > {1/C)xn' for all n. Finally, the terms ln(x^'°/x^'=) on the right hand 
side of (|4.7p are uniformly bounded in n and k. Combining these ideas shows that 
for n large enough 



:^ki{tn)xl" 



In 



„yki 



y,yki 



Cfel 



p 

EE 

i = l {i— 



In 



^Vk 



< 



-qi.nX' 



Vkl 



(4.10) 



where gi,„ — > cxd as n — ^ oo. Thus, we have used one half of a single term enumerated 
in (|4.8I) to bound every term enumerated on the right hand side of (|4.7p . 

Since we already know that the terms in (|4.8|) are all strictly negative for large 
enough n, all that remains to show is that the terms in (14. 9p . which are all positive 
for large enough n, are also dominated, in a similar manner as (j4.10p . by some terms 
in (|4.8I) . as n ^> oo. 

Pick a reaction from (|4.9I) . yo ~^ Voj say. Suppose that the source of the reaction 
is a complex in tier T^, and the product is in tier Ti-m for some m > 0. By the weak 
reversibility of the network, there is a series of reactions beginning with j/g and ending 
with j/o such that no reaction is enumerated more than once (that is, there is a path 
along the directed diagram with no yk y'f. used multiple times). We now claim 
that there must be a subset of these reactions, enumerated as ri,. . . ,rb, satisfying 
the following conditions (below, we denote the tier of the source complex for reaction 
as g and the tier for the corresponding product complex as T^j p): 

1. di_s < i — m. 

2. For £ 6 {1, . . . ,b}, the source complex of is in a strictly higher tier than 
the corresponding product complex; that is, di^s < di^p. 

3. For ^ > 2, we have di^s £ di-i,p- 

4. db^p > i. 

5. For £ e {1, . . . , 6}, we have di^s < i- 

Note that, for example, the series of reactions ri, . . . ,rb above can be constructed 
from the original series by only taking the first b reactions for which the source is in a 
strictly higher tier than the product, but stop (i.e. pick b) once a product complex is 
in a tier that is equal to or lower than that of j/q. If yo — ^ J/o is a reversible reaction, 
then we may take 6=1 with the reaction ri simply being the reverse reaction 2/q — > j/o 
from tier Ti_m to tier Ti. 

By condition 2 we know that each of the reactions ri , . . . , rt, are enumerated in 
(|4.8|) . Let do = max^jd^^s}, and let ydg G Td„ be a choice of complex from tier Td„. 
Note that do < i hy condition 5, and so 



Xn''° /x"^" -> OO, as ri OO. 



(4.11) 



By construction and an application of the triangle inequality, there is a constant 
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Ci > such that for 71 large enough (so that the In terms dominate the c^^ terms) 



In 



> rjCiXn' 



In 



1=1 



(4.12) 



where the apparent "flip" of the ratio conies from taking the absolute value of a 
negative term, and where 77 > is the parameter used to bound all the functions 
Kk{t): r] < Kk(t) < for all t > 0. By condition 2 in the construction above we 
have that each term in the product on the right hand side of (|4.12p goes to +00, as 
n 00, and hence the entire product does. Further, by conditions 1, 3, and 4 above, 
there is a C2 > such that 



nX-n 
y 

1=1 x„ 



> G 



r.yo 



uniformly in n. We now may conclude that there is a C3 > so that 
■ / ^ y^e\ b 



i=l X„ 



>C3\n{x'f;>/xl«) + co, 



(4.13) 



for all n. Let M be the number of reactions enumerated in ()4.9p over alH S {1, . . . , P}. 
Combining (|4.13p with (|4.12p and (|4.1ip . we see that for n large enough 



2M '^''^ {tn)Xn ' 

e=i 



In 



+ Cr 



In 



^ya 
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(4.14) 

where (72. n 00, as n ^> cxd. Note that the first term on the left hand side of (|4.14p 
is the left hand side of (|4.12p divided by 2M, and the second term is the summand 
of (|4.9p associated with the reaction yo — > yg. As yo — > j/g was arbitrary, we may 
conclude that for n large enough. 
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< 



which arc the terms of (|4.8p plus the terms of (|4.9p . 

Combining the inequalities (|4.10l) and (|4.15p shows that for n large enough, the 
summation found on the left hand side of (|4.7I) . or equivalently on the left hand side 
of ()4.6p . must be strictly negative. This is a contradiction with ()4.6p holding for all 
71. Therefore, we must have that condition C2 holds. □ 

Lemma 4.8. Let {S,C,TZ^K,{t)}, with S =^ {Si, . . . , Sn}, be a non- autonomous 
mass-action system with bounded kinetics and a single linkage class. Suppose xq G 
M^Q is such that a;o) remains bounded and dist((/)(i, xq), 9M>q) — as t — )■ 00. 
Then, there does not exist a subsequence of times t„ — 00 such that C is partitioned 
along Xn — (f>{tn,Xo) in which Ti consists of a union of linkage classes. 

Proof. Note that in the one linkage class case Ti can only consist of a union of 
linkage classes if Ti = C. We suppose there is such a sequence of times, tn — > 00, such 
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that C is partitioned along a;„ (t>{tn, xq) with Ti = C (and there are no other tiers). 
By the boundedness of the sequence, we may consider a convergent subsequence with 
hmit point z. Let U = U{z) = {z S {1, . . . , : zt = 0}. Note that U is nonempty 
because the trajectory goes to the boundary. Note that C is necessarily partitioned 
along this subsequence as well, with the same tiers. By Theorem 14.61 there is a non- 
negative conservation relation w G ]R>q that respects the pair (J7, {Ti}). 

Because the support of w is C/ ^ 0, and (j)i{tn,Xo) — >■ for all i S ?7, we have 
that w ■ (t){tn,xo) ^ 0, as n — )■ cxd. However, because Ti = C, we also have that 
'u^ ' {y'k~ Vk) = for all Uk ^ y'k ^ T^- Thus, we see from (|2.2|) that w ■ (j){t, xq) > is 
constant in time t, contradicting that w ■ xq) — > as n — >■ oo. □ 

Lemma l¥!5] savs that condition C2 of Lemma l4.7l can never hold in our setting, and 
so CI must. The following, final lemma will allow us to conclude that any trajectory 
satisfying such a family of Lyapunov functions must converge to a single point. 

Lemma 4.9. Let {S,C,TZ,IC{t)}, with S = {Si, . . . , Sn}, be a non- autonomous 
system with bounded mass-action kinetics. Suppose xq G is such that for any 
X G ]R>o, there exists a T = T^ > such that t > T implies 

l^^(^W) < 0, 

where x{t) — (j){t,X()) is the solution to the system with x{0) — xq and kinetics IC{t). 
Then a;(0(-,a;o)) is a single point. 

Proof. Note that the trajectory remains bounded because each Vx{x{t)) does. 
Also, we have that for any x G M>q there exists a e^- > such that 

Vx{x{t)) — > C-, as i ^> oo, 

where the non-negativeness of Cx follows by the fact that Vx{x) > for all x G K>o- 
The boundedness of x{t) implies there is at least one w-limit point of the trajectory. 
The question now is: can there be more than one? Suppose so. That is, we assume 
the existence of 21,22 G '^(4>i'jXo)) with zi ^ Z2. 

Note that for any x G R>o, Vx{zi) — Vx{z2) = Cx, where if Zi G 9IR->o we define 
Vx on the boundary via its continuous extension to the boundary. Let xi,X2 G M^q 
be arbitrary. Then, after some algebra we have 

= VxAzi) - VxAz2) - {VxA^.l) ~ Vx,{z2)) 

= [zi - Z2) ■ (ln(x2) - In(xi)) . 

But, Xi,X2 G K>o, and hence (lna;2 — Inxi) G M^, were arbitrary. Thus, zi = Z2. □ 
We now have our main result. 

Theorem 4.10. Let {S,C,TZ,IC}, with S — {Si, . . . , S\^\}, be a weakly reversible, 
single linkage class chemical reaction network with mass-action kinetics. We assume 
that for xo G M!>o trajectory (j){t,xo) satisfies the following two conditions 

1. (j){t,Xf)) is bounded (int), and 

\s\ 

2. u}{(j){-,XQ)) is either completely contained in dWy^ or completely contained 
within the interior of M^fg . 

Then a;((/)(-, xq)) H 9M>q = 0, and the trajectory is persistent. 

Remark 1. Note that the conclusion of the theorem guarantees that uj{(j){- ,xq)) 
is completely contained within the interior of the strictly positive orthant. Said difjer- 
ently, uj{(I){-,xq)) can not be contained within dWy^. 
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Proof. Suppose, in order to find a contradiction, that for this xq there is at least 
one z e a;(0(-, xo)) fl 9M>q. Let 

U = {i & {1, . . . , \S\} : Zi — for some z e oj{(f>{-, a^o))}: 

which is nonempty as there is at least one z e 9R>q. That is, these are ah the indices 
for the species whose concentrations approach zero along some subsequence of times 
for this specific, fixed trajectory. Therefore, and equivalently, i e C/ if and only if 

\imiid (j)i{t, Xq) — and linisup (/>i(t, xq) < cxd, 

where the second fact follows from the boundedness of trajectories, whereas for j ^ U 
we have 

< liminf xo) < limsup0j(t,a;o) < oo- (4.16) 

Let {Su,Cu,Ti-u} denote the reduced reaction network of {S,C,TZ} associated 
with U (see Definition 13. II) . and let IC{t) — ICu{t) denote the projected dynamics (see 
Section [3.2p . with Kk(t) denoting the non- autonomous variables defined via (13. 7p that 
take the place of the rate constants in standard mass-action kinetics. It is important 
to note that by (|4.16p and the definition of the Kfc(t)'s given in p.7p . we have the 
existence of an 77 > such that 

r] < Kk{t) < 1/77, (4.17) 

for all t > and ah k e {1, . . . , |7^|c/|}. That is, {Su,Cu,Uu,K:{t)} is a non- 
autonomous system with bounded mass- action kinetics. Another way to see this 
fact is just to note that the chemical species whose concentrations are uniformly 
bounded from above and below, those j ^ U, have been incorporated into the Hkit), 
thereby yielding a bound like (|4.17l) . By Lemmas 13.31 and 13.41 the reduced network 
{Sij,Cij,TZij} is weakly reversible and has only one linkage class. 

We let \Su\ — N and denote by x{t) e M^q the solution to the reduced dynamical 
system for this specific trajectory. By condition 2., above, which pertains to the 
original system, and by the construction of U, the set of oj-limit points of the trajectory 
of the reduced system must exist on 9R>q . Combining Lemmas 13. 3[ 14. 7[ 14. 8[ and 14.91 
shows that the set of tj-limit points of the trajectory of the reduced system, x{t), 
must consist of a single point. By construction of U, this point must be the origin 
£ M^, as otherwise there is an i £ [/ for which liminft_^oo Xi{t) > 0, a contradiction 
with the definition of U. However, we also know by Lemmas 14.71 14.81 and 14.91 that 
-^Vx{x{t)) < for t large enough, where X £ R^o is 

arbitrary. Therefore, because the 
origin is a local maximum of Vx, we can not have that x{t) <E M^. □ 

The following corollary, which was the main goal of the paper, states that the 
Global Attractor Conjecture holds in the single linkage class case. 

Corollary 4.11. Let {S,C,TZ,JC} denote a complex-balanced system with one 
linkage class. Then, any complex-balanced equilibrium contained in the interior of a 
positive compatibility class is a global attractor of the interior of that positive class. 

Proof. Trajectories of complex-balanced systems satisfy conditions 1 and 2 in the 
statement of Theorem 14. 10) [TJ. The result then follows by the discussion in Section 
11.11 after the statement of the Global Attractor Conjecture. □ 

In particular, if {iS,C,7^, /C} is a weakly reversible, deficiency zero system with a 
single linkage class, then the conclusion of the Global Attractor Conjecture holds. 
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Note that the single hnkage class assumption in Theorem 14.101 and hence Corol- 
lary mTTJ was only used in conjunction with Theorem 14.61 in the proof of Lemma 
to guarantee that the top tier, Ti, could not consist of a union of linkage classes. If 
it can be guaranteed in any other way that the top tier, in the construction detailed 
in the previous lemmas, can not consist of a union of linkage classes, then the conclu- 
sions of Theorem 14.101 and Corollarv l4.1H that complex-balanced equilibria are global 
attractors of their positive classes, will still hold. Also, note that if it can be shown 
that condition 2 of Theorem 14.101 is always satisfied by weakly reversible networks 
with mass-action kinetics, something we believe to be true, then the Persistence Con- 
jecture, as stated in Section [01 of this paper, will also be proven in the single linkage 
class case by the arguments in this paper. 

For completeness, we state the following corollary to the proof of Theorem 14.101 
which points out that the rate constants are permitted to be bounded functions of 
time before the projection. 

Corollary 4.12. Let {S,C,K,IC{t)}, with S = {5*1, . . . , ^i^i}, be a weakly 
reversible, single linkage class chemical reaction network with bounded mass-action 
kinetics. We assume that for xq G R!>q the trajectory (f>{t,xo) satisfies the following 
two conditions: 

1. (j){t^XQ) is bounded (int), and 

\s\ 

2. a;(0(-,xo)) is either completely contained in 9M>q or completely contained 

within the interior of R^q . 

Then a;(0(-,xo)) H 9K>q = 0, and the trajectory is persistent. 

Proof. The proof is exactly the same as the proof of Theorem 14.101 □ 

5. Example. We present one example to demonstrate the ease with which the 
main results can be applied. Consider the system with reaction network 



which was recently used as an example of a system whose persistence was beyond the 
scope of known theory 125] . This network has four complexes, one linkage class, and 
the dimension of the stoichiometric subspace can be checked to be three. Thus, the 
deficiency is zero. As the network is clearly weakly reversible, the Deficiency Zero 
Theorem implies that, regardless of the choice of k^, the system is complex-balanced. 
Theorem l4.10l and Corollarv l4. 1 11 tell us that the system is persistent, and that each of 
the complex-balanced equilibria are global attractors of their positive stoichiometric 
compatibility classes. Thus, the long term behavior of the system is now completely 
known. 

Acknowledgements. Thanks to Masanori Koyama for a thorough reading of an 
early draft that found an error in the proof of Lemma 14.71 Also thanks to two 
incredibly meticulous reviewers who helped with the presentation substantially. 
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